Odd spin-triplet superconductivity in a multilayered superconductor-ferromagnet 

Josephson junction. 
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We study the dc Josephson effect in a diffusive multilayered SF'FF'S structure, where S is a 
superconductor and F,F' are different ferromagnets. We assume that the exchange energies in the 
F' and F layers are different {h and H, respectively) and the middle F layer consists of two layers 
with parallel or antiparallel magnetization vectors M. The M vectors in the left and right F' layers 
are generally not coUinear to those in the F layer. In the limit of a weak proximity effect we use 
a linearized Usadel equation. Solving this equation, we calculate the Josephson critical current 
for arbitrary temperatures, arbitrary thicknesses of the F' and F layers {Lh and Lh) in the case 
of parallel and antiparallel M orientations in the F layer. The part of the critical current IcSR 
formed by the short-range (SRC) singlet and S = triplet condensate components decays on a 
short length = \/D/H, whereas the part IcLR due to the long-range triplet IS] = 1 component 
(LRTC) decreases with increasing Lh on the length = sjDj-nT . Our results are in a qualitative 
agreement with the experimenlj^. 

PACS numbers: 74.45.-l-c, 74.50.-|-r, 75.70.Cn, 74.20.Rp 



I. INTRODUCTION 

According to the Bardeen, Cooper and Schrieffer— the- 
ory of superconductivity in conventional metals and their 
alloys the superconducting condensate consists of singlet 
Cooper pairs. These pairs can be described by a wave 
function / which is in the absence of the condensate flow 
symmetric in the momentum space (s-wave singlet pair- 
ing) . In terms of the creation and annihilation operators 
ip'^^tp, this function can be represented in the form of 
a thermodynamic average fsngit — t') ^ (^t(0''^i(^') ~ 
V'i(^')V't(0)- At equal times {t = t') this function de- 
termines the order parameter A : A = A(V'-f(OV''i(^)) — 
— X{ipl{t)ijj-f{t)) , where A is the coupling constant of an at- 
tractive interaction. The spin-independent scattering by 
ordinary impurities does not affect this type of supercon- 
ductivity (Anderson-theorem)^. 

In the last two decades other types of superconductiv- 
ity have been discovered. The most important example 
is high Tc superconductivity discovered by Bednorz and 
Miiller-2. in cuprates. 

As a result of intensive study of these superconductors, 
it was demonstrated that, although the Cooper pairs in 
this case are also singlet, their wave function essentially 
depends on the momentum p and changes sign with vary- 
ing momentum direction in CuO planes. In the simplest 
version, the dependence of the order parameter A on p 
has the form: A(p) = Ao{cos^ {pxa^) — cos"^ {pyUy)) {d- 
wave singlet pairing). Such a dependence A(p) allows 
one to construct the so called tt— Josephson junction con- 
sisting of single crystals of high Tc superconductors with 
an appropriate orientation of these crystals with respect 
to each otheii^. The ground state of this junction cor- 
responds to the phase difference equal to tt. 



Another type of superconductivity (triplet) has been 
discovered in strontium ruthinate Sr2RuCu4^^'^ and in 
heavy fermion intermetallic compounds*. In contrast to 
the singlet superconductivity, the wave function of the 
Cooper pairs /f-|-(p -I- p\t — t') ~ (?/'-|-(p,t)V't(p'ji')) is an 
odd function of momentum p, so that for equal times 
t — t' the function /^-|-(p,0) and the order parameter 
A(p) do not equal to zero (p-wave triplet pairing). Only 
for some directions of the momentum p the order param- 
eter turns to zero. This means that, in agreement with 
the Pauli principle, the pair wave function changes sign 
under permutation of spins and momenta. The momen- 
tum dependence of the condensate function makes the 
singlet d-wave and triplet p-wave superconductivity sen- 
sitive to scattering even by potential (not acting on spin) 
impurities. 

An unusual mechanism of superfluidity was proposed 
by BerezinskyS,. Having in mind liquid He^, he consid- 
ered a retarded interaction between atoms and assumed 
that the order parameter A(w) and the wave functions 
/^-|-(p,aj) or /^_|,(p,a;) in the Matsubara representation 
were even functions of momentum but odd function of 
the Matsubara frequency uj. However, experiments on su- 
perfluid He^ revealed that the p-wave triplet type of the 
superfluidity is realized in He^ rather than the one pro- 
posed by Berezinsky^^'^^ . Some possibilities to realize the 
exotic Berezinsky-type mechanism of superconductivity 
in various systems in context of the pairing mechanism 
in high Tc superconductors were considered in Kefsr^^'^. 

This exotic type of superconductivity (or superfluid- 
ity) was regarded for quite a long time as a hypothetical 
one. Only recently it has been realized^'^ that the odd 
triplet superconductivity might exist in a simple SF bi- 
layer system consisting of a conventional s-wave singlet 
BCS-type superconductor S and a ferromagnetic layer F 
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with a nonhomogeneous magnetization M. 

To that moment it had already been well known that 
in an SF system with a homogeneous magnetization M 
the Cooper pairs penetrated the ferromagnet over a short 
length ^i? — D/E(,x (in the diffusive limit), where 
D = vl/3 is the diffusion coefficient, I = vt is the mean 
free path, and E^x is the exchange energy in the fer- 
romagnet (we set the Planck constant h equal to 1). 
Since the exchange energy usually is much larger than 
the critical temperature of the superconductor Tc, the 
length is much shorter than the length of the conden- 
sate penetration into a normal metal in an SN bilayer: 

The Cooper pairs penetrating the ferromagnet with an 
uniform magnetization consist of electrons with opposite 
spins. Their wave function is however the sum of a sin- 
glet and triplet component with zero total spin projection 
on the z-axis {S = 0). The exchange field mixes these 
components and the triplet component with S" = is 
unavoidable in the ferromagnet. The sum of these two 
components can be considered as a short-range compo- 
nent (SRC). 

The part corresponding the triplet 5* = component 
has the form Unit " t') ~ {Mt)Mt') + 
At equal times t — t' this function equals zero in agree- 
ment with the Pauli principle. Therefore the function 
ftr-fiit ~ t') is an odd function of the time difference 
(t — t') or Lo in the Matsubara representation. The order 
parameter in the superconductor S is related only to the 
singlet function /sng(w) which is an even function of w. 
The superconducting order parameter in F is zero if the 
coupling constant Ai? in the Cooper channel equals zero. 

The situation changes if the magnetization orientation 
in the vicinity of the SF interface is not fixed. This case 
was analyzed in Ref.^^, where an SF bilayer with a do- 
main wall located at the SF interface was considered. It 
was shown that in such a system not only the singlet and 
triplet S = Q components but also the odd triplet compo- 
nent with iS* = ±1 arises in the ferromagnet. The latter 
component penetrates the superconductor over a large 
distance that does not depend on the exchange field and 
is of the order ^^r provided the spin- dependent scattering 
is not too strong. 

This odd triplet component can be considered as the 
long-range triplet component (LRTC). As the LRTC is 
symmetric in the momentum space, the scattering by po- 
tential impurities does not affect this component. 

In subsequent theoretical papers various types of SF 
structures where the LRTC may arise were studied (see 
review articlesi^"— and references therein) . In Ref the 
creation of the LRTC is predicted in a diffusive SF struc- 
ture with a Bloch-type domain wall (DW) . The width of 
the DW, w, was assumed to be larger than the mean free 
path I : I <^ {w, ^f}- 

A more general case of the DW with a width, arbi- 
trary with respect to the mean free path, in a SF struc- 
ture with an arbitrary impurity concentration was stud- 
ied in Ref.2i. The LRTC in diffusive SF structures with 



a Neel-type DWs has been analyzed in Ref.--. The case 
of a half-metallic ferromagnet in SF or SFS structures 
was investigated in Refsj^i^i^. Braude and Nazaro\*2^ 
studied the LRTC in SF structures with a highly trans- 
parent SF interfaces so that the amplitude of the conden- 
sate functions induced in the ferromagnet was not small 
(strong proximity effect). Ballistic SF structures with a 
nonhomogeneous magnetization, where the LRTC could 
be created, were studied in Refs.^'*"^''. The papers^^"— 
were devoted to the study of the LRTC in spiral ferro- 
magnets attached to superconductors. 

In several paper o^^'^"^!^^" — the LRTC was investigated 
in SF structures with the so-called spin-active interfaces. 
In the approach used in these papers, the properties of 
the SF interface are characterized by a scattering matrix 
with elements considered as phenomenological parame- 
ters. In this approach one does not need knowing the 
detailed structure of the SF interface and can proceed 
calculating physical quantities using these parameters. 
We wiU see that even in the framework of the quasiclas- 
sical theory one can obtain effective boundary conditions 
for the LRTC provided the width of the DW w attached 
at the SF interface is thin enough (such an approach was 
used in Refj^). From the physical point of view the re- 
gion with a narrow DW can be regarded as a spin-active 
SF interface. If the width w is comparable with the Fermi 
wave length, one has to go beyond the quasiclassical the- 
ory and derive the boundary conditions from the first 
principles (see^ as well as^i^ and references therein) . 

By now, several papers presenting a quite convincing 
experimental evidence in favor of the existence of the 
LRTC have been published22i^i26i37 j^gf 32 ^-^e con- 
ductance of a spiral ferromagnet {Ho) attached to two su- 
perconductors was measured. It was concluded that the 
conductance variation below the superconducting criti- 
cal temperature is too large to be explained in terms 
of the singlet component. Keizer et al.'^^ observed the 
Josephson effect in an SFS junction with a half-metallic 
ferromagnet Cr02- The thickness of the F layer was 
much larger (up to ^1 mkm) than the penetration depth 
of the short range condensate components. Moreover, in 
the metal where free electrons with only one spin direc- 
tion are allowed, no pairs with opposite spins are possible. 
Therefore, only triplet l^l = 1 component can survive in 
this ferromagneti^i. However, there is no controllable pa- 
rameter in this system that would allow one to change the 
amplitude of the LRTC. A similar long-range Josephson 
effect in a SFS junction with Cr02 as a ferromagnetic 
layer was observed by Anwar et al. in a recent work'^''. 

Recently the dc long-range Josepshon effect has been 
observed in a more complicated SFS structure with a 
controllable parameter'^*'. In the experimental setup of 
this work, F was not a single ferromagnetic layer but a 
multilayered structure of the NF'NFNF'N type, where N 
is a nonmagnetic metal, F' is a weak ferromagnet {PdNi 
or CuNi) and F is a strong ferromagnet (Co). The mid- 
dle F layer was in its turn a trilayer structure consisting 
of two F layers with antiparallel orientation of magneti- 
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zation M and of a thin layer (Ru) that provides RKKY 
coupHng between the F layers. 

The authors of Ref.— measured the Josephson critical 
current Ic for different thicknesses L of the F' and F layers 
(we denote the thicknesses of the F' and F layers as 
and Lh layers respectively). It was demonstrated that in 
the absence of the F' layers {Lh = 0) the critical current 
Ic was negligible if the width of the F laye r Lh essentially 
exceeded the small length S^h — \/L)/H, where H is the 
exchange energy in the F layer. This is what one expects 
for the conventional superconductivity. However, adding 
the F' layers resulted in an increase of the critical current 
Ic by several orders. The dependence of Ic{Lh) is non- 
monotonous: the critical current is small at small and 
large Lh reaching a maximum at Lh ~ ^/i- 

The authors of Ref.^^ suggested an explanation of these 
results in terms of the LRTC. Note that the mean free 
path / in the structure studied in Ref.— is rather short 
(the diffusive limit in the F' layers and an intermediate 
case in the F layer). 

Theoretically the dc Josephson effect in multilayered 
SFS junctions with a non-collinear magnetization orien- 
tation has been studied in several works. In Refsi^i^S 
the Josephson current in ballistic SFS junctions was 
calculated. The diffusive SFF'S junctions were consid- 
ered in Refs.— However the long-range Josephson 
effect in the junctions with two F layers is not possi- 
ble; the Josephson critical current Ic is not exponen- 
tially small only if the total thickness of the ferromagnetic 
layer, Lp + Lpi , is comparable with the short length ^p: 
Ic ^ exp(-(LF + Lf,)/£,f). 

The diffusive Josephson junctions with three ferromag- 
netic layers and non-collinear M orientation, where the 
long-range Josephson coupling may exist, have been ana- 
lyzed in Refs42.i^. The authors of RefJ^ considered the 
F'SFSF'structure with different magnetization M orien- 
tations in the F and F' layers. In Ref.™ a somewhat 
different, but more suitable for experimental realization, 
SF'FF'S structure with different M directions in the F 
and F' parts was analyzed. In both papers the exchange 
energy in the F and F' was assumed to be equal. 

The amplitude of the LRTC, /i, and the Josephson 
critical current due to this component IcLB. are calcu- 
lated in both works. Although the structures studied in 
Refs.— !^ are different, the results obtained are similar. 
The final formula for the critical current can be written 
in both the cases as 

IcLR. = F{Lh) sin ai sin ar (1.1) 

In Eq. (jl.ip . ai^r are angles between the z-axis and 
the magnetization vectors in the left (right) F' layers, 
while the magnetization M in the F layer is assumed to 
be parallel to the z-axis. The function F{Lh) is a non- 
monotonic function with a maximum at Lh ^ ^h- This 
function vanishes at small and large thickness Lh of the 
F' layers (see Eq. (12) in Ref.'^^ and Fig. 2 in RefM). 

Qualitatively, this prediction agrees with the observa- 
tions in Ref.'^^. However, experimental parameters pre- 



sented in this publication are well defined and this makes 
a more detailed comparison of the theoretical predictions 
for LTRC with the experimental results quite interesting. 

In this paper we analyze an SF'FF'S structure which, 
being formally similar to that considered by Houzet and 
Buzdin^, is in many respects different and closer to the 
structure studied experimentally. 

First, unlike Refj^, we assume that the exchange en- 
ergies in the F' and F layers {h an H respectively) are 
different. 

Secondly, we assume that the SF' interface is not per- 
fect and the proximity effect is weak. This assumption 
allows us to linearize the Usadel equation and to calculate 
the critical Josephson current Ic at any temperatures T 
(in Ref.^'^ only the case of temperatures close to Tc was 
considered). 

Thirdly, we also analyze the case when the F layer con- 
sists of two domains with parallel and antiparallel orien- 
tations of magnetization. 

At last, we derive a formula for the current Ic for arbi- 
trary thicknesses, Lh^H, of the F' and F layers (in Ref.''^ a 
formula for Ic is presented only in the limit Lh/^h ^ !)■ 

The paper is organized as follows. In Section |ll] we 
formulate the problem and write down necessary equa- 
tions. In Section IhII we analyze the case of thin F' layers 
{Lh/^h ^ 1)- We calculate the amplitudes of the short 
range singlet and triplet S — components as well as the 
LRTC. In Section Hv] the case of arbitrary lengths Lh.H 
will be considered under assumption that the angles ai^r 
are small. Using the formulas obtained for the ampli- 
tudes of different components, we calculate in Section Ivl 
the critical current Ic in the limiting cases: a) for par- 
allel and antiparallel M orientations in the F layer and 
for arbitrary angles ai^r assuming that the thickness Lh 
is small, b) for arbitrary Lh^n under assumption of small 
ai^r- The results obtained are discussed in Section IVll 



II. MODEL AND BASIC EQUATIONS 

We consider a multi- layer S/F Josephson junction 
shown schematically in Fig.l. It consists of two supercon- 
ductors, S, and three ferromagnetic layers F, F'^j. The 
middle F layer may consist of two domains or layers with 
parallel or antiparallel orientations of the magnetization 
M. A similar Josephson junction has been studied ex- 
perimentally in a recent work'^^. 

The presence of normal N layers in the experimental 
S/NFiNFNFrN/S structures cannot change qualitatively 
the results for the S/F'iFF'i/S structure obtained here 
because the scattering in the N layers does not depend 
on spin (if a weak spin-orbit scattering can be neglected) . 
Therefore, all the superconducting components, singlet 
and triplet, decay in the N layers in a similar way over 
a large distance of the order ^jy. The exchange fields 
acting on electron spins are h in the F' layers and H 
in the middle F layer. The magnetization vector M in 
F is supposed to align along the z-axis and it has the 
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FIG. 1: (color online) Josephson structure under considera- 
tion. The F'(F) layers are weak (strong) ferromagnets. The 
middle F layer consists of two layers with parallel or antipar- 
allel (shown in figure) magnetization orientation. The arrows 
show the directions of the magnetization in F' and F layers. 



components M (0, sina/^^, cos ai^r) in the F\r layers. The 
magnetization in the F layer is oriented along the z-axis 
but may have parallel or antiparallel orientations in the 
regions {—Lh < x <Q) and (0 < a; < Lh)- 

For explicit calculations we use the quasiclassical 
Green's function technique, which is the most efficient 
tool for studying SF structures (see reviewsii"— i^"— ) , 
and assume that all the ferromagnetic layers are in the 
diffusive regime, so that the Usadel equation can be used. 
The amplitude of the condensate wave function in the 
ferromagnetic layers is assumed to be small (weak prox- 
imity effect) and therefore the Usadel equation can be 
linearized. The smallness of the condensate wave func- 
tion is either due to a mismatch of the Fermi velocities 
in S and F or due to the presence of a tunnel barrier at 
the S/F interfaces. 

The anomalous (Gor'kov) quasiclassical Green's func- 
tion in the considered case of a spin-dependent interac- 
tion is a 4 X 4 matrix /. We are interested in the dc 
Josephson current /c, i.e. in a thermodynamical quan- 
tity. Therefore we can use the Matsubara representation 
for the matrix / and consider / as a function of the 
Matsubara frequency uj = 7rT(2n -I- 1) and coordinate x 
normal to interfaces: / = f{uj,x). The linearized Usadel 
equation for / has the form (see^, Eq.(3.15)) 



d'^jjdx^ - kU - i{Kj,/2) cosa(x) 
x{tana(x)f3«)[(T2,/] + [0-3,/] + } ^ 



(2.1] 



where = 2\bj\/D, k\ = 2sgnuj ■ h/D in the F'l^r layers 
and — 2sgnuj ■ H/D in the F layer, the Pauli ma- 
trices Ti^di operate in the particle- hole and spin space 
respectively. The angle a is equal to ai^r in the F'l^^ lay- 
ers and to zero in the F layer in the case of the parallel 



orientation of the magnetization M in the domains. In 
the case of the antiparallel orientation a{x) = tt in the 
interval {—Lh < x < 0) and a{x) — in the interval 
(0 < X < Lh)- The diffusion coefficient D is assumed to 
be the same in all the ferromagnetic layers. 

The matrix / can be represented for the system un- 
der consideration in a form of an expansion in the spin 
matrices tr, as 



(2.2) 



The matrices (Tq and ai^s are the unit matrix and the 

ax,z Pauli matrices, respectively. The /o,i,3 matrices are 
matrices in the particle-hole space. The first term is the 
short range triplet component with the zero projection 
of the total spin on the z-axis, the second term is the 
LRTC with the non-zero projection of the total spin, and 
the third term is the singlet component of the condensate 
Green's function (seei^ii^). 

Eq. (|2.ip should be complemented by boundary con- 
ditions. We consider the simplest model of the S/F het- 
erostructures assuming that the interfaces have no effect 
on spins (spin-passive interface). These boundary condi- 
tions have the forr" ^^'^^ 



(2.3) 



df/dx\x=±L = ±lBfs\x=±L, 



where jb = ^/{Rb<^), Rb is the S/F interface resistance 
per unit area, a is the conductivity of the ferromagnet. 
The matrix fs is the Gor'kov's quasiclassical Green's 
function in the left and right superconductors. It has 
the form 



fs\x=±L = /s<5'3'8'('r2COS(/3±fiSini^), 



(2.4) 



where fs — A/y/uj^ + A^ , ±1^9 is the phase in the right 
(left) superconductor, so that the phase difference is 2ip. 

If there is a spin-dependent interaction in a thin layer 
at the interface (exchange field, spin-dependent scat- 
tering etc), the boundary condition acquires a more 
complicated form. In particular, the coefficient 75 be- 
comes a matrix with matrix elements containing very 
often unknown phenomenological parameters. Such in- 
terfaces are called spin-active interfaces. In many pa- 
pers the LRTC is studied in SF systems with spin-active 
interfacesS^^i^^. 

The F/F'i^r interfaces are assumed to be ideal and 
therefore the function f{x) and its derivative df/dx must 
be continuous at these interfaces. Solving the linear equa- 
tion (|2.ip with the boundary conditions (|2.3I) one can cal- 
culate the dc Josephson current using the formul a^^'^^ 



= i<j/8)2TTTY,Tr{a„®fJdf/dx}, 



(2.5) 



j>0 
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This problem can be solved in a general case of an 
arbitrary thicknesses of the F and F' layers {Lu and L^) 
and angles ai^r- 

However, the general results are too cumbersome. In 
order to present analytical formulas in a more or less 
compact form, we consider two limiting cases: a) thin 
F\i- layers {Lh ^ £.h,£.N) and arbitrary angles a^,; , b) 
arbitrary thicknesses L//, Lf^ but small angles a^,; (a <C 
1). In the next section we consider the case a). 



a; > 0, (3.6) 

f±{^) — ^± cosh(K;j-a;) cos iy9 + Bj- sinh(K;j-a;) sin iy9, 
a; < 0, (3.7) 

The relations between coefficients A, A and B, B 
should be found from the continuity of the matrices /± {x) 
and their derivatives df±ldx at a; = 0. This gives: 
A — A, khB — RhB. Using the boundary conditions 
(|3.4l3.5p , we find the coefficients A± , B± 



III. THIN F' LAYERS 

In this section we assume that the F'l r layers (or h- 
layers) are very thin so that the inequality \Kh\Lh ^ 1 is 
satisfied, where k\ = 2sgnuj(h/ D) (usually kh ^ i^n = 
\J TxT jD and therefore the condition unLh ^ 1 is also 
fulfilled) . In the zero order approximation in the param- 
eter KhLfi, the exchange field in the entire ferromagnetic 
region except the thin F'l ^ layers is homogeneous and 
equal to H. Thus, only the /o,3 components in the ex- 
pansion (|2.2p that describe the short-range components 
(SRC) differ from zero. 

The matrix / satisfies the equation 

d^f/dx^ - nlf - in^j cos a{x)ai®f = 0, (3.1) 

where k\j = 2sgnuj{H/ D). 

The angle a{x) =0 in the case of parallel orientation 
of the vector M, and a{x) = at x > 0, a{x) = tt at 
a; < in the case of antiparallel orientation. We rewrite 
Eq. p.ip for the diagonal in spin space components f± — 
fo ± /s as 

d^f±/dx^ - 4/± =0, a: > 0, (3.2) 
d^f±/dx^ - kJ± =0, x<0, (3.3) 

where ± the case of the parallel 

M orientation in both domains (a: > and x < 0) and 
— if magnetization vector at x < 

changes sign with respect to its direction at x > 0. 
The boundary conditions for matrices f± follow from Eq. 

df±ldx = ±-1 B fsih cos ip + Ti sirup), 

X - iff, (3.4) 

df± /dx = Tlsfsih cos (fi - Ti sin (fi), 

X = -Lh; (3.5) 

The solution of Eas. (|3.2j|3.3l) can be sought in the form 
f±{x) = A-t cosh(K±x) cos iy9 -|- i?± sinh(K±x) sin iy9. 



A± = ±^[fs{LH)cosh0± + fs(-LH)cosh0±], 
B± = ±2±[fs^LH)—sinhe± - fs{-LH)sinhe±], 

V± K± 

(3.8) 

where 2?± = k± sinh 6± cosh 9± + R.± sinh 9± cosh 9± , 
0± = k±Lh and 0± = k±Lh- 

In the case of parallel (P) and antiparallel (AP) orien- 
tations of the magnetization in the F layer we obtain the 
function 'D± 

T>±p = 2k± sinh 6'± cosh 6*^ , 
Vj^ap = X'-AP = 2i?e(K+ sinh0-|_ cosh6'_). 

(3.9) 

One can see from Eqs. (|3.6ll3.8p that the SRC de- 
cays exponentially away from the SF interfaces over the 
short length £^h- Indeed, at \Lh — x| ^ £,h we obtain 
from Eqs. (|3.6ll3.8p that /±(x) ^ {jb / K±)fs{T2 cosip + 
fi siiiip)exp{-{LH - x)/^h)- 

Let us now find the LRTC. First, we obtain the ef- 
fective boundary conditions for this component. Assum- 
ing that Lh <^ £,h, we can integrate Eq. (|2.ip over the 
thickness of the F'-layers and come to effective boundary 
conditions for the triplet component 

dFi/dx\x^±LH = ±71/3 (±Lh) sin Qfr,;, (3.10) 

where 71 = nf^Lh- 

We have introduced in Eq. p.lOp a matrix Fi = T^®fi 
describing the LRTC. This matrix Fi satisfies an equa- 
tion that directly follows from Ea. (|2.1l) 

d'^Fi/dx^ - kIFi = 0. (3.11) 
The solution for the matrix Fi can be written as 

Fi = Ai cosh(K(jx) -f- Bi sinh(Ktjx) (3-12) 
From the effective boundary conditions (I3.10p we find 
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"^1 = n [/a (-^ff) sin a,. + /a (-iff) sin a/], 

-^1 = t; — — [/3(iff)sinar - /3(-iH)sina;]. 

(3.13) 

where 71 = 2sgnuj{h/ D)Lh,, i.e. the matrix Fi is an odd 
function of the Matsubara frequency. 

The solution for Fi, Eq. p.l2p . demonstrates that the 
LRTC described by the function Fi decays slowly at a 
large distance of the order ^ ^jv- The matrix /a in 
Eq. is expressed through /± : /a = (/+ - /-)/2, 

where the matrices f±{±LH) are given by Eqs. p.6ll3.8p . 

As it should be, the function Fi{x) turns to zero in 
the absence the exchange field or in the case of collinear 
magnetization because 71 = Kf^Lh ^ h and sina^.i = 
in the case of collinear M orientations. 

Eqs. p. 6113. 8|3. 12113. 13]) determine all the condensate 
Green's functions. Using these functions, we calculate 
the Josephson current in Section fVl 

IV. ARBITRARY THICKNESSES OF 
FERROMAGNETIC LAYERS AT WEAK 
NON-COLLINEARITY 

Consider now a more interesting case of an arbitrary 
thicknesses of the ferromagnetic layers F', F (or h, H - 
layers). We restrict ourselves with the case of the parallel 
M orientations in the F layer because there is no quali- 
tative difference between the behaviour of the LRTC in 
the P and AP magnetic configurations. For simplicity 
we assume that the angle a is small, a <C 1. In this 
case the amplitude of the LRTC is proportional to the 
small parameter a. In the zero order approximation only 
the singlet component, /a , and the short range triplet 
component, /o, with zero projection of the total spin of 
Cooper pairs on the z-axis are not zero. Indeed, we will 
look for a solution of Eq. (12. ip in the form /3(a;) ^ 
foix) fi{^) ^ {cosh(Ka;), sinh(Kx)}, where n is the 
eigenvalue. 

In the ferromagnetic layers we obtain the following 
equations for the eigenvectors 

fo{n^ ~ kI) ~ hin^cosa = (4.1) 
/3(k^ — K^) — /qik^ cosa — sina = (4.2) 
Fi{k^ - + sin a = (4.3) 

where the matrix Fi introduced in Eg. p.lOp describes 
the LRTC. This set of equations has three eigenvalues 

2 2 2i-2 2 2 I A A\ 

'^1,2 = I^F± = '*cu ± ^/^Fj '^a = i^u, (4-4) 

Two of them, kf± , describe a sharp decay of the 
density of Cooper pairs in the ferromagnet (in the case 



{H,h} ^ T,A) and the latter one, k^^ {k^j = l/^Lfl), is 
an inverse characteristic length of decay of the LRTC in 
the ferromagnet. By order of magnitude it is equal to 
« TiT /D, which shows that the length is rather 
large and does not depend on the exchange energies ft., H. 
Spin-orbit interaction or a spin-dependent impurity scat- 
tering make this length shorter— li^iS"— 

4 = 1^1+ «m (4-5) 

where « Ta.\\i{DTra,DTsp-orb}, Tm and Tsp-orb are 
characteristic times related to the spin-dependent im- 
purity scattering or spin-orbit interaction. The lengths 
Kp]j_ also depend on and can be found by shifting 

It is seen from Eqs. (|4.m4.3p that the LRTC arises only 
at non-zero a when Fi ^ Q. In the zero-order approxi- 
mation {a = 0) we should find the matrices /o.a in each 
ferromagnetic layer. As follows from Eqs. (|4.m4.2p . at 
a — only the eigenvectors /o± = ±/3± corresponding 
to the eigenvalues k± can be finite. 

The solution for /o,3 satisfying the boundary condi- 
tions (|2.3p can be written as 

f'ii^) = y^AT2AHi COshjKHjX) COS (fi 
j 

+fiBHj smh{KHjx) sin ip} (4.6) 

fo{x) = y^(-l)J+i{f2^gj cosh(Kgjx) cosip 
j 

+TiBHj sinh{KHjx) sin ip} (4.7) 
in the i?-region (F layer) and 

h{x) = ^{HChj cosh{KhjS:) + Sl^^ smh{Khjx)] 
j 

+Ti [Cg^ cosh(K,„i) + smh{Kh,S:)]} (4.8) 

fo{x) = J2i~^y^'i^^^(^S cosh(^,,i) + Sl^^ sinh{nf.,5:)] 
j 

+fi [Cj^^ cosh(K,y i) + s'h] smh{Kh,S:)] } (4.9) 

in the /i-region (F' layers), where a: = a; — L, j = l,2so 
that khi,2 = kh± etc. 

The coefficients AHi^Bui and C^h\ , s'^^\ , c''^] , S^^\ are 
found from the boundary conditions (|2.3p . We write 
down here the expressions for Cui , Sui (see Appendix) 

Ah± = \±s'il/R±{cs-sc), (4.10) 
Bh± = X±sl2/R±icc;ss) (4.11) 

where X± = [Kh/KH)± and R±{cs\sc) = 
(cosh sinh 61/^ + Asinhfl/i cosh6'/f )±, R±(cc]ss) = 
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{cosh 9 h cosh 0H + Asinh^/i sinh0^f)±. The coefficients 



^(1,2) 



are equal to 



(1) 

h± 



iBfs 



sin if\ 5 



(2) 
h± 



isfs 

2Kh- 



■ cos if 



(4.12) 



Eqs. (14.6114. 12p determine tfie SRC in tliree-layer 
Joseplison junction with, different exchange energies in 
the middle (//-region) and terminal (/i-region) F layers. 
We use Eas. (|4.6ll4.7p for the calculation of the Josephson 
current due to the SRC. 

Let us turn to the calculation of the LRTC. We write 
the equation for the matrix Fi in the /i-region projecting 
of Eq. (j2.ip on the ai matrix in the spin space 

d^Fiix)/dx'' - KlFiix) = ~4sma ■ fsix), (4.13) 



where the function /3(x) in R.H.S is given by Ea. (|4.8l) . 

The solution of Eq. (|4.13p can readily be obtained (see 
Appendix). In the //-region the function Fi{x) obeys 
the same equation but without the R.H.S. {a — 0). The 
solution in this region has the form of Eq. p.l2p . 



V. JOSEPHSON CURRENT 

Using the Green's functions, /o,i.3 obtained in Sections 
mil and IIVI one can now calculate the dc Josephson cur- 
rent. Substituting the expansion (|2.2p into Eq. (|2.5p . we 
obtain for the Josephson current density 



JJ 



iaiTT V Tr{T3[fodfo/dx+f3df3/dx+f,dh/dx]}, 



j>0 



The first two terms in Eq. ()5.1|) are the contribution 
from the SRC (/ca) and the third term is due to the 
LRTC (/i = T3F1). The part of the Josephson current 
density which is caused by the SRC, jcSRi can be written 
also in the form 



3.JSR = I'^^T Tr{h[Udf+/dx + Udf^/dx]}, 

uj>0 

(5.2) 

where f± — fa i fs are the diagonal matrix elements in 
spin space. The part of the Josephson current density 
which is caused by the LRCT, jcLR, can be presented as 



Fi = AiH cosh^Ki^x) + BiH smh{Ki^x). 



(4.14) 



In order to calculate the Josephson current Ij we need 
to know the coefRcients Aim and Bin- Considering the 
symmetric case, ar — ai — a, we obtain (see Appendix) 

AiH = f2AiHC0Sip, BiH = fiBiHSinLp (4.15) 
where 



Ai/sinh^^ =isina ^ (-l)^+'^g^ 

J = l,2 

s'^hj, sinhe Hj . , „ XjCoshOHj. 

x{^ cosh6';ia;^-7 ^ - 1 + smh ^/^^ ^— f }; 

Uj[cs;sc) Kj(cs;sc) 



BiHCOshe^ ^isma ^ (-1)^+^5'^^^ 



(4.16) 



3 = 1.2 



X { ^ [cosh ehu> - 1] + sinh Ohu; \ ^^^^ }; 

Kj[cc;ss) Rj{cs\sc) 

(4.17) 

where 9^ = k^L, 9,-,^ = Ku^Lh, 9hj = i<-HjLh (ki,2 = 
K±). The functions Rj{cs]sc) and Rj(cc]ss) are defined 
in Eqs. (j4. 10114. lip . 

In the limit of thin /i-layers (|0/i±| <C 1), we see that 
the products Am sinhfl^j and Bin cosh^^j agree with the 
coefficients Ai and Bi in Eq. p.lSp . This means that 
the amplitude Fi of the LRTC goes to zero at Lh — > 0. 
On the other hand, it is seen from Eos. (|4. 16114. 17p that 
at \9h±\ 3> 1, the functions Rj are exponentially large 
and therefore the amplitude of the LRTC decreases with 
increasing the thickness Lh- 



J 



JLR = -iaTrTY,Tr{f3[FidFi/dx]}, (5.3) 



j>0 



We use these formulas for the calculation of the critical 
Josephson current in the limiting cases a) and b). 

A. Thin F' Layers 

Consider the limit of thin F' layers. Substituting 
Eqs. ([3?6l) into Eq.(|01). we obtain 



3JSR = Tr{f3[K+A+B+ + (5.4) 

CJ>0 

With the help of Eq. (|3.8l) this expression can written 
in the form 



jjSR = JcSfl sin(2(^), 

jcSR 



(5.5) 



= a7|7rTV/J[-^ + -^], (5.6) 



ui>0 + 



Using Eas. p.9p . we find the critical current density due 
to the SRC for the P and AP magnetic configuration 



jcSR,AP = cr7|27rT Y fi 



smh{29H+) 
1 



2/?e(K+ sinh 9h+ cosh 9h-) ' 



(5., 
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FIG. 2: (color online) Normalized critical current due to the 
SR components, IcSR, for the parallel (P) and antiparallel 
(AP) orientations of magnetization in domains in the mid- 
dle _ff-layer vs the thickness of the iif-layer (A) and tem- 
perature (B). The thickness of the /i-layers is assumed to be 
small (Lfi << ^ft). The normalized temperature T, exchange 
energy H and thickness Lh are measured in units Ao and 
= y^D/Ao, where Ao = A(r) at T = 0. The lower and 
upper curves correspond to the P and AP configuration, re- 
spectively. The critical current jc is measured in units of its 
value in a junction with _ff = (no exchange field) and nor- 
malized length Lh = 0.5. The point and upper thin curves 
correspond to H = 70; the lowest and upper thick curves cor- 
respond to H — 30. The normalized temperature (in Fig. A) 
and thickness Lh (in Fig.B) are equal to 0.1 and 0.5; the 
parameter r^A equals 1. The values of the critical current 
corresponding to -ff = 30 are reduced by 10 times. 

These formulas show that the critical current den- 
sity jcSR decays exponentially with increasing the length 
Lh over a short scale of the order ^h- The de- 
crease of the current density jcSR.P with increasing 
the thickness Lh or exchange energy H is accompa- 
nied by oscillations^^"^^. Oscillations in the dependence 
jcSR.Ap{LH, H) are absent in the case of antiparallel ori- 
entations. The latter behavior was predicted by Blanter 
and Hekking^^ who considered antiparallel orientation of 
magnetization in ferromagnetic domains in a SFS Joseph- 
son junction and presented formula for jcSR.,AP in the 
limit 9h+ > 1. Eq.dlH) generahzes Eq.(26) of Ref.^^ ^j. 
the case of arbitrary 9h±i i-e. arbitrary Lh- A rapid 
decay and oscillations of the critical Josephson current 



FIG. 3: (color online) The same dependence as in Fig. 2 for 
the normalized critical current due the LRTC components, 
IcLR- As in Fig. 2, the thickness of the ft-layers is assumed 
to be small (L^ << ^h). The normalization units are the 
same as in Fig. 2. The curves are numbered so that: 1 - AP, 
{H = 30); 2 - AP (70); 3 - P (30); 4 - P (70). 



in SFS or SIFS junctions were observed in many experi- 
mental works^i^"— . This oscillatory behaviour was pre- 
dicted a few decades ago^^iZi. 

In Fig. 2a and 2b we show the dependence of 
the normalized critical currents density IcSR = 
jcSR{h, H, LH)/jciO,0, Lh) originating from the short- 
range component (the components /o,3) on the thickness 
of the i/-layer Lh and temperature for the P and AP 
magnetization orientation. The plots are obtained on 
the basis of Eqs. (|5.7ll5.8p applicable in the case of thin 
/i- layers. We see that in the case of the P orientation the 
critical current density jc,SR caused by the SRC changes 
sign with increasing Lh, while in the case of the AP ori- 
entation the current jc,SR is always positive. However, at 
a fixed Lh the current density jcSR decays monotonously 
with increasing temperature. For a smaller exchange en- 
ergy the dependence jcSR{T) has another form and may 
change the sign (see the next subsection). 

Note that the critical current for the AP orientation 
IcSR,AP is always larger that for the P orientation. The 
critical current LcSr.ap in an SFIFS Josephson junction 
with the antiparallel magnetization orientation in the F 
layers may be even exceed the critical current in SIS junc- 
tion without the F layers provided that the coupling be- 
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FIG. 4: (color online) The dependence of the normalized criti- 
cal current due to the SR component (point line) , LRTC (solid 
thick line) and the total critical current (solid thin line)on the 
thickness of the Jf-layer for H = 70 (A) and H = 170 (B). 
The normalized temperature is equal to 0.1. Other parame- 
ters are the same as in Fig. 2. 



tween S and F layers is strong enough ( here I stands for 
an insulator) 

The Josephson current Jjlr due to the LRTC is found 
using Eqs. ([01 13.12|3.13l) . We obtain 

. ^ 2 . . V- Tr{Tsh{LH)f3hLH)} 



w>0 



Kuj smh{29 Hoj) 



(5.9) 

The matrices f^liLn) are expressed in terms of 
the matrices f±{±LH) : fsi^LH) = [/+(±Lh) - 
/+(±Ljj)]/2. Using Eqs. p.6ll3.7p . we find that jjlb. = 
JcLfl sin(2(^). The critical current jcLR iu the case of the 
P and AP orientations is given by 



jcLR,p = -4cr7B7i sinarSinQ!;(27rr) 



E 



/I 



K„ sinh(20^f;^) 



RePi ■ ReP2, 



(5.10) 



^>o '^\p'^^ sinh(26'//„) 



[Pi - 1] 



(5.11) 



jcLR.AP = -criBli sina^ sina/(27rT) 



where Pi = cosh^ 61+ /!?+ , P2 = sinh^ 6*+ /!?+ , P3 = 
|cosh6'+|2(l + |tanh6'+pi?e(K+/K_)), V+ and Vap = 
Vap+ = Vap- are defined in Eq. dSJ]). Eq. (|530)) 
resembles Eqs. (12) and (10) of RefsJ^ii^, respectively. 
Comparing Eq. (|5.10p with Eq.(lO) of Ref4^ one has to 
have in mind that other boundary conditions and tem- 
peratures close to Tc are considered in that work. Note 
that the parameter (ji/kuj) = iKhLh){Kh/ Kuj) may be ar- 
bitrary; we assumed that KhLh <C 1, but the parameter 
(Kh/n^) is large. 

In Figs. 3a and 3b we plot the same dependencies as in 
Figs. 2a and 2b for the critical current jcLR. originating 
from the LRTC. This current decays monotonously with 
increasing both the temperature and Lh- Note that the 
critical current jcLR has the opposite sign to jcSR.AP if 
sin ar sin a; > and in the case of the AP configuration 
is always larger than in the case of the P configuration. 

In Figs. 4a and 4b we plot the most important de- 
pendencies of the jcSR{T) and jcLR{T) as well as the 
total critical current, jcTot = jcSRiT)+ \ jcLRiT)\, on the 
length Lh for different magnitudes of the exchange en- 
ergy H. The parameter p{ar, ai) = 47^^^ sina^ sin a/ is 
taken to be equal to 0.2. The value oi H = 170Ao ap- 
proximately corresponds to the exchange energy in Co 
(i?ea; = 309mey, see^*) used in the experiment^. One 
can see that for H — 70 Aq and H ~ 170Ao the critical 
current is caused by the LRTC at Lh ^ O.S^a and at 
Lh > 0.4^A, respectively, where = v^-D/Aq. The lat- 
ter curve is close to the one observed in Refi^ provided 
we accept ^a = lOnm corresponding to the value of the 
diffusion coefficient D about IQcm? / s. 

One can say that the obtained results are in a qualita- 
tive agreement with experimental datap^. It is difficult to 
carry our the quantitatvie comparison because our model 
is simplified. We admit the standard model in which one 
neglects the difference in the diffusion coefficients for the 
major and minor electrons in the ferromagnetic layers al- 
though this difference in Co is large. We also assume the 
diffusive limit for all layers. At the same time, the mean 
free path in the strong ferromagnet (Co) seems to be 
larger than . In this case the formula for the Josephson 
current even for one-layer SFS junction becomes rather 
complicated^^. 

In addition, the conductivities in all the layers are as- 
sumed to be equal whereas in experiment the conduc- 
tivities in layers (Cu, Co and PdNi) are different. The 
interface resistance strictly speaking, is not known 
and can be only estimated. However, one can see from 
Eqs. I|5. 7115.5)) that both the "effective conductivity" a 
(averaged over all layers) and the averaged interface re- 
sistance Rb enter the expression for the critical current 
density jc as a prefactor before the sum. It disappears 
when we plot the critical current normalized to its value 
ai h = H = Q. The most interesting dependence of the 
critical current jc on temperature and thicknesses L^^h 
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is determined by exponential functions. 

We assumed that the proximity effect is weak, that 
is, the amphtude of the condensate functions in fer- 
romagnets is small. As follows from Eq. (|3.8l) . this 
means that the parameter |7b/k+| \pLh,H / Rb(^+\ ~ 
\pLh.H/ Rb\ ^ 1- According to Ref^S, where a struc- 
ture similar to that in Ref.'^*' (but without strong ferro- 
magnets) was studied, the interface resistance per unit 
area, Rb = SRsf, was equal Rb — 2.3 • lO^^ficm^ 
whereas RpS = ppLp ~ 7- 10~^^f2cm^, that is, the ratio 
RfS/Rb « 10""* is indeed very small, here S is the cross 
section area of the junction. Taking into account that the 
resistance of the h and H {Co) - layers are comparable, 
we conclude that the proximity effect in experiments^ is 
weak. 

As to the value of the critical current, we do not 
attempt to carry out a quantitative comparison with 
experimental value because it depends on the ratio 



{pL/Ri 



{RpS / Rb) which is known only on the 



order of magnitude. In addition to that, the Joseph- 
son junction used in experiment contains many interfaces 
each of which reduces the proximity effect. 



B. Arbitrary Thicknesses of the F Layers 

In order to calculate the current density jjsR, we sub- 
stitute Eas. (|4.6ll4.12p into Eg. dO)) using the relations 

fo,3{^) — [f+{^) ^ /+(2;)]/2. After simple calculations, 
we obtain for the critical current density jcSB. in the case 
of the P configuration 



jcSR = a(2^r)7| V Re[ 



/I 



^ ^''kh+R{cs/sc)+R{cc/ss)+ 



(5.12) 

The critical current density due to the LRTC, jcLR, for 
the case of the arbitrary thicknesses Lh.H is found from 
Eas. (|4. 14114. 1715.5]) . For the symmetric case (a^ = o;; = 
a ^ 1) we get 



cLR = -a{2TTT)-^%a^ ^ 



/I 



sinh(20tj) 



where 



Im{M)Im{N) 
(5.13) 



M 



sinh 9h+ cosh Kuj cosh 9h+ sinh 

R{cs/sc)+ kh+R{cs/ sc)^ 



(5.14) 



^ _ cosh 9h+ cosh 6*/,^^ ^ sinh 6*^+ sinh dhu 



R{cc/ss)^ 



kh+R{cc/ ss). 



(5.15) 



and 6cj = 0ijh + ^uH- In the antisymmetric case (a^ = 
— a; = a ^ 1) the sign of the critical current density 



should be changed. Thus, in this case the critical current 
has the same sign as in SNS junction. 

Using Eas. (|5.5ll5.15p . we calculate numerically the crit- 
ical currents jcSR and jcLR for the case of symmetric 
case {ar = ar = a <^ V) and the P configuration. 
In Fig. 5a and 5b the dependence of jcSR. and jcLR 
on temperature T is shown for the case when the ex- 
change energies are not very high {H — 2h = IOAq) 
and the lengths Ln.h are small {Lh ~ Lh — O.I^a)- It 
is seen that the current jcSR depends on T in a non- 
monotonic way and can change sign. This type of a 
non-monotonic temperature dependence was obtained in 
many works, both experimentai^ and theoretical (see 
review o^^'^^ and references therein as well as a recent 
paper—). The critical current due to the LRTC jcLR 
decays with increasing T monotonously. Note that a 
non-monotonic behavior of jcLR{T) was found in Refi^ 
for a half- metallic ferromagnet. In our model, we do 
not find values of parameters at which this dependence 
would not be monotonic. However, there is no contra- 
diction between these two results because, as was shown 
recently^, the non-monotonic temperature dependence 
of the critical Josephson current IcLR takes place only 
for a sufficiently large exchange energy (comparable with 
the Fermi energy ep)- In our study we assume that both 
exchange energies, h, H, are small in comparison with ep- 
In a recent work^^, the Josephson dc effect was observed 
in a SFS junction with ferromagnetic Cu2MnAl-]ieus\eT 
alloy as a ferromagnetic layer. In the interval of thick- 
nesses 8nm < Lp < 10.5nm the critical current as a func- 
tion of Lp shows a very slow decay with increasing Lp, 
and its temperature dependence was a non-monotonic in 
this interval of Lp. 

In Fig.6A we show the dependence oi jcLR{Lh) for dif- 
ferent exchange energies h setting the angle a equal to 
VO.l. It is seen that the critical current caused by the 
LRTC has a maximum at Lh ~ S,h, that is, the maximum 
shifts to smaller with increasing h. The same depen- 
dence jcLR{Lh), but for different H, is shown in Fig.6B. 
One can see that the position of the maximum weakly 
depends on H. A nonmonotonic behaviour of the critical 
current as a function of Lh was observed in experiments^. 



VI. DISCUSSION 

We have considered the long-range triplet component 
in an SF'FF'S diffusive Josephson junction with a non- 
coUinear magnetization AI orientation in the F' and F 
ferromagnetic layers. Assuming that the proximity effect 
is weak, we have solved the linearized Usadel equation 
and found the pair wave functions for the short-range 
(singlet and 5 = triplet) and long-range components 
(the LRTC with \S\ — 1) in the cases of different ex- 
change fields in F' and F layers, arbitrary temperatures 
and parallel (antiparallel) M orientation in two domains 
in the middle F layer. 

Our study was motivated by recent experimental re- 
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FIG. 5: (color online) Temperature dependence of the normal- 
ized critical current due to the SR (A) and LRTC (B) com- 
ponents for the P orientation. The upper and lower curves 
correspond to h = 5, H = 10 and to normalized lengths 
Lh ~ 0.1 and Lh — 0.12. The point curve corresponds to 
h^5,H = 10 and Ln = 0.1. The "length" Lh = 0.1 and the 
parameter TmAo = 0.1. 



suits concerning the observation of the LRTC in a multi- 
layered (seven layers between superconductors) Joseph- 
son junctionsSS. The model used in our study, although 
accounting for the properties of the SFS junction used in 
the experiment^, is somewhat simplified. In particular, 
we solve the Usadel equation in which the difference in 
transport properties of the minor and major electrons in 
the F layers was ignored. This approximation is rather 
crude especially for strong ferromagnets (in Co the dif- 
fusion coefficients for the minor and major electrons may 
differ by an order of magnitude^^). Account for differ- 
ent transport properties for electrons with up and down 
spins leads to a change not only in the boundary condi- 
tion p.Sl! ^'^ but also in the collision integral by potential 
impurities'"^. The scattering by stationary fiuctuations 
of magnetic moments in the ferromagnetic layers is also 
taken into account in the simplest way ignoring the spin- 
orbit interaction^. 

Thus, the approximations taken by us do not allow a 
quantitative comparison of the critical current jc mea- 
sured in experiment and calculated on the basis of the 
quasiclassical theory within a simplified model. The 
formulas for jc, Eas. (|5.7ll5.8|) . contain an effective con- 
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FIG. 6: (color online) A) Normalized critical current due to 
the LRTC as a function of the thickness of the /i-layer for 
different exchange energies h for the P orientation. The pa- 
rameters are: h = 2,5, 20 (solid thin, solid thick and point 
lines, respectively). The other parameters are H = 50, T = 
0.25, TmA = 0.1, Lh = 0.5. B) The same dependence for 
different H: H = 5 (thin solid curve), H = 20 (thick solid 
curve) and H = 100 (point curve). The parameter h = 5. 
Other parameters are the same as in Fig.6A. All the curves 
correspond to the P orientation. 



ductivity averaged over all layers and effective interface 
resistance. Even in a simpler case of SFS junctions with 
a single ferromagnetic layer it is not possible yet to ob- 
tain a quantitative agreement between theory based on 
the quasiclassical theory and experimental^. 

Fortunately, both these parameters enter the corre- 
sponding formulas as prefactors which disappear when 
the critical current Ic{T, h, H, Lh, Lh) is normalized to 
its value a,t h = LI — 0. Therefore, the performed 
study allows one to understand what kind of depen- 
dencies of the critical current Lc on different parameters 
(T, Lh, Lh, LI, h) can be obtained in multilayered Joseph- 
son SFS junction. If the magnetization in the domains 
in the F layer is aligned parallel, the critical current IcSR 
caused by the short-range components oscillates with in- 
creasing the thicknesses of the F' and F layers [Lh^n)- 
This behavior was predicted in Refs.— and observed 
in many experimental works^i^"— i^"— . The formulas 
obtained in this paper generalize previous theoretical re- 
sults for IcSR (see reviews^— and references therein) to 
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the case of different exchange fields in the F' and F films. 

If the magnetization in domains in the F layer are an- 
tiparallel, the critical current IcSR decays exponentially 
in agreement with the results of Rcf.^^. In both the cases 
of the parallel and antiparallcl orientations the character- 
istic length for the decay of IcSB. is of the order of £,h,H- 

The critical current IcLR due to the LRTC decays in 
both cases with increasing Lh and T in a monotonic way 
on the length of the order of ^at. For certain values of 
parameters of the system the total critical current Ic = 
IcSR + IcLR coincides with IcSR ^f small thickness of the 
F layer Lh and with IcLR at larger Lh (see Fig. 4). This 
behavior agrees with the dependence Ic{Lh) observed in 
the experiment^. 

As it was predicted in Ref."'^ for a F'SFSF' system and 
in Ref.'*'^ for a SF'FF'S system, the current IcLR has a 
maximum as a function of at Lh ^ decaying to 
zero at large and small Lh- Thus, the presence of the 
the F' layers makes the F layer with a strong ferromag- 
net "transparent" for the LRTC. Unlike previous theo- 
retical studies, we present a formula for the amplitude 
of the current IcLR for arbitrary thicknesses Lh^H, tem- 
peratures and exchange fields h,H. This allows one to 
choose optimal parameters for observing the LRTC. 

However, it remains unclear how the angle dependence 
of the critical jcLR{oir,i) shows up in the critical current Ic 
observed in experiment. Generally speaking, orientations 
of the M vectors in F' and F films are not necessarily the 
same {ar.i 7^ 0). If there are domains in the F' layers, the 
total critical current Ic,Av is determined by averaging the 
current density jc over the width of the junction as the 
magnitude and sign of the average (sin sin ai) depends 
on the number of domains, orientations of the magneti- 
zation M in these domains etc. This issue requires a 
further theoretical and experimental studies. 



VII. APPENDIX 



Substituting Eqs. (|4.9l4.8p into Eas. ((7T|) . we obtain 

Eqs. (|4T2)) . The coefficients C^j;' are found from the 

matching conditions of the functions fo.aix) and their 
derivatives ed x = Lh- They equal 



(1) _ oil) R±{SC/CS) ^ ^(2) _ ^(2) R±{SS/CC) 



h± 



h± 



R± {cs/sc) 



(7.2) 



Eqs. (|4. 8114. 12117. 2|) determine the form of the function 
/3(x) in Eq. (|4T3l) . The solution for Eq. (|4T3l) is given 
by the formula 



Fi{x) = Fiunix) + Fi„„„(a;), 



(7.3) 



where i^i„„ (x) , Finun (x) are solutions of a homogeneous 
Eq. ()4.13|) (without the R.H.S.) and Finunix) is a partic- 
ular solution of this equation. These functions are 



Fiunix) = T2[a*-^-' cosh{Ki^x) + 6^^-* smh{Ki^x)] 
+fi [a^^^ cosh(K^i) + 6^^) sinh(K^i)]} (7.4) 



Rj{cc/s 



Finunix) = sm a {■^2<S;j^'[-— — - — -cosh[Khj 

JrL'i [ CS I SC] 

R ' ( SC I cs\ 

-I- sinh(K/yx)] -I- — - — - cosh(K;yx) -I- sinh(K;ija;)]} 



(7. 



Matching Eqs. (17. 3117. 5p . we find finally the matrices an 
and 6h, Eqs. (|4.16M.17|) . 



The boundary conditions for the matrices /0.3 follow 
directly from Eq. (|2.3p and have the form 
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